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SUMMARY 


A multiple- layer radiative transfer model of a vertically inhomogeneous, 
optically shallow water mass is briefly described. This model is directed 
toward use in remote sensing of water properties. Some preliminary results 
and qualitative predictions are presented. 


INTRODUCTION 


In most applications of remote sensing involving water quality the 
assumption is made that the water is vertically homogeneous. Usually the 
water is also assumed to be optically deep, i.e. absorption and reflection 
by the bottom are taken to be negligible. These assumptions are frequently 
adequate, as evidenced by the wide-ranging success in using remote sensing 
for observation of water properties. However, some concern has been voiced 
with respect to the general validity of the standard assumptions 
(ref. 1) and, in at least one case, a changing vertical distribution of 
material in water has been linked to variation in remote observations 
(ref. 2). It is the underlying thesis of the work presented here 
that vertically inhomogeneous and/or optically shallow waters are fairly 
common, that the inhomogeneity will affect the remotely sensed upwelling 
radiance and, therefore, that there is need for a mathematical model appli- 
cable to these situations and useful for remote sensing applications. 

In the following pages a multiple-layer radiative transfer model of an 
optically shallow water mass is briefly described. In order for this model 
to be directly useful in remote sensing applications it must be invertible. 
This requirement necessitates several simplifying assumptions which will 
inevitably limit the accuracy of the model in at least some situations. 
Hopefully the advantages to be gained by having an easily manipulated model 
of a rather complex system should outweigh the loss in accuracy. Initially 
it is intended only that the model give a good qualitative description of the 
system, although care has been taken to formulate the model in such a way as 
to facilitate using the model to make quantitative predictions. 
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DESCRIPTION OF THE MODEL 


Assumptions and Definitions 

Radiative transfer models tend to be extraordinarily complex and utterly 
resistant to inversion. Much of the complexity arises because of the diffi- 
culty in describing anisotropic light fields. While the anisotropy is an 
important part of the interaction of light with water, it need not be described 
in minute detail. Characterization of the light field in general terms should 
provide the simplification necessary for deriving a model capable of des- 
cribing a vertically inhomogeneous system but still susceptible to inversion. 

The basic assumption is that the underwater light field may be effec- 
tively characterized by the apparent optical properties: the diffuse atten- 

uation coefficients for upwelling and downwelling irradiance, the irradiance 
reflectance, and the radiance reflectance. The diffuse attenuation coefficient 
for downwelling irradiance is defined as 


-1 dE^(z) 
" E (z) dz 

Q 


( 1 ) 


where Ej(z) is the downwelling irradiance at depth z. Likewise, the diffuse 
attenuation coefficient for upwelling irradiance is defined as 


-1 E^(z) 

k (z) = — — p — X -z 

u Ej(z) dz 

a 


( 2 ) 


The diffuse attenuation coefficient will be dependent to some extent on 
the radiance distribution. However, evidence is growing that the underwater 
radiance distribution does not vary in a way that strongly affects the diffuse 
attenuation. In fact, observations by Baker and Smith (ref. 3) indicate that 
the diffuse attenuation coefficient for downwelling irradiance is remarkably 
insensitive to the radiance distribution, whether due to changes in sun angle 
or depth. 

The irradiance reflectance, R(z), is defined as 


R(z) 


E^(z) 

Ed(z) 


(3) 
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The term R(z) , which has become a fairly standard measure of water color 
is seemingly independent of illumination conditions (refs. 4 and 5). We 
may also define a radiance reflectance: 




L^(z,e’) 


(4) 


where L^(z,0*) is the upwelling radiance at depth z and in direction 0’. 


In addition to the above standard properties we define two other para- 
meters: the irradiance scattering function for downwelling light, B^(z), and 

a single-scattering irradiance attenuation coefficient for upwelling light, 
k^(z). defined as the irradiance dE (z) scattered upward at depth z 

from a horizontal slab of thickness dz when iYluminated by the downwelling 

irradiance, E - (z) : 
cl 


dE (z) 

" E (z)dz 
cl 


(5) 


The scattered irradiance, E (z), is attenuated as it proceeds toward the 
water surface. However, the irradiance attenuation coefficient for upwelling 
light, k^(z), is not appropriate since it implicitly includes the backseat tar- 
ing of downwelling light already described by B^(z). To adjust for this we 
define a single-scattering irradiance attenuation coefficient, k' (z) , such 
that: 

k’ (z) = k (z) - B (z) (6) 

u d 

Like the other apparent properties, B^(z) and k' (z) will be assumed to be 
quasi-inherent optical properties sincie they are dependent on the radiance 
distribution in essentially the same way as k^, IL , and R. 


Model Geometry and Final Equations 

This model treats the water as a plane-parallel medium of arbitrary depth 
in which the optical properties, depth and thickness of each layer, as well as 
the depth of the water and the bottom reflectance, may be specified independ- 
ently. Figure 1 illustrates the model geometry for the relatively simple 
situation of one layer of turbid water in an otherwise homogeneous water 
column. The attenuation coefficients of pure water are k and k^, 
while the irradiance scattering function for pure water is B’ . The correspond- 
ing optical properties of homogeneously distributed substances in the water 

are n k , n k’ and n B' where n is a concentration parameter which may vary 
ss^ssss s 
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between 0 and 1. For n = 0 the concentration of the material is zero. For 

g 

n = 1 the concentration of the material is a maximum, i.e., the water has 
s 

been replaced entirely by the substance. 

An intermediate layer is shown in figure 1. This layer contains material 
which is optically distinct from the surrounding water. The optical properties 
of the material present in the layer are complete 

analogy to the homogeneously distributed material. 

A portion of the irradiance above the water surface, E^, is reflected at 
the air water interface, (S • The remainder is transmitted ([1-S^]E^). 

The irradiance is attenuatet as it passes down through the water column. At 
depth z, a portion of the downwelling irradiance is scattered back toward the 
water surface. This scattered irradiance is further attenuated as it travels 
to the surface and across the air-water interface. 

A portion of the downwelling irradiance reaches the bottom (z = d) and 
is diffusely reflected there (^^) • The reflected portion is attenuated as it 
travels upward through the water column. The irradiance is affected by the 
local optical properties at each depth which are assumed to be constant 
within each layer. 

The irradiance reflectance immediately below the water surface is given 
by 


-2(K.-K )Az. r -K (h.+Az.) -K h. 


R(O-) = 


B -2(K.-K )Az. r 


o J 1 
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where: k = (1-n ) (k + k^ ) + ti (k + k') 

o s w w s s s 

B^ = B^(z) for homogeneous regions 

m = number of intermediate layers 


h. 

1 

Az, 

1 

K. 

3 


depth to the top of layer i 
thickness of layer i 


B. 

3 

k. = 


kl 

X 


(1-n -n.) (k +k’) + n (k +k' ) + n. (k.H-kO 
s X ww s s S X X X 

Bj(z) for h.<z<h.+Az 
d XX 

downwelling irradiance attenuation coefficient for layer i 
upwelling irradiance attenuation coefficient for layer i 
concentration parameter for layer i 


and where Az^ = and first term on the right hand side of 

equation (7) describes the portion of R^(0') due to scattering from all the 
areas in which only the homogeneously dxstributed material is present; the 
second term describes the return from each of the layers; the third term 
accounts for the bottom reflectance. In deriving equation (7) it was assumed 
that none of the layers overlapped. 


The irradiance reflectance is related to the observations immediately 
above the water surface by: 


(0) - p (6) L (0) Q(0’) n , 

U cl K, W 

(1 - S )E E * 1-p (0’) 

a d r w 


where : 


L^(6) = upwelling radiance above the water surface in direction 0 
L^(0) = downwelling sky irradiance in direction 0 
p (0) = specular reflection at the water surface in air 

3 . 

p (0') = specular reflection at the water surface 
w 

S = diffuse reflectance of the air-water surface 

a 

n = index of refraction of water 
w 

= downwelling irradiance above the water surface 

= portion of the upwelling irradiance internally reflected in 
the water 
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Q(0^) = conversion factor relating radiance and irradiance reflec- 

tance 


The above equations were derived in analogy to single-scattering radiative 
transfer models and as such will be called a single-scattering irradiance 
model (SSI). The primary distinction between the (SSI) model and two flow 
theory is that the downward scattering of upwelling light is ignored. This 
simplification causes results for strongly scattering waters to be inaccurate. 
In spite of its appearance, equation (7) is a relatively crude representation 
of the reflectance characteristics of an optically shallow, vertically inhomo- 
geneous water body. The simplifying assumptions used in deriving equation (7) 
will limit the absolute accuracy of the model; however, it should provide a 
good qualitative description of variations in ocean color. Moreover, this 
model may be accurate enough in some situations to yield moderately accurate 
quantitative predictions. 


APPLICATION OF THE RADIANCE MODEL TO AN 
OPTICALLY SHALLOW, HOMOGENEOUSLY ABSORBING WATER COLUMN 

Upon formulation of the radiative transfer model for the case of an 
optically shallow, homogeneously absorbing water column, measurements 
of volume reflectance were conducted under the controlled conditions of a 
water tank. In an attempt to maintain a simple and inexpensive experimental 
design, sunlight was utilized as the illumination source. Variations in the 
absorptive capacity of the water as well as in column depth were considered. 


Model Formulation 

In applying the radiance model to the case of an optically shallow, 
homogeneously absorbing water column, a number of simplifications may be 
applied to equation (7). Each simplification is based upon one or more of 
the following assumptions: 

1) the water column contains no intermediate layers possessing 

unique optical properties (Az = 0) 

2) the concentration of any particulate scattering material sus- 

pended throughout the water column is very small 

3) the bottom is highly reflective and closely resembles a 

completely diffuse reflector 

4) internally reflected irradiance, E^, is negligible 

5) observations are made in the nadir direction 0=0 
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Combining equations (7) and (8) and applying the above assumptions results in 
the relationship 


L 

u 


Experimental Apparatus 

Experiments were conducted using the cubic meter tank illustrated in 
figure 2. The sides of the tank were painted with a low-reflectance, ultra- 
flat black paint; the floor of the tank was coated with a high reflectance, 
flat white paint in order to optimize the return signal. Shadows from the 
sides of the tank were avoided by conducting all the measurements at a time 
when the sun was highest. At no time was the solar zenith angle greater than 
35°. 


(1-p^) (1-S ) 


'w' a' 4. T 

■ + Pk\ 


Q n 2 
w 


(9) 


Radiance measurements were made using a United Detector Technology 
Spectral Radiometer designed to continuously scan the visible and near 
infrared portion of the light spectrum from 400nm to llOOnm. As shown in 
figure 2, the radiometer was positioned directly over the center of the tank 
so as to record upwelling radiance in the zenith direction. In this config- 
uration the radiometer shaded the water surface directly underneath from 
downwelling sky radiance. Therefore, the reflectance term the 

right side of equation (9) may be neglected. 

The total downwelling solar irradiance was measured with the use of a 
panel coated with a standard reflectance medium, barium sulphate. Such a 
coating is noted for its high reflectance and close resemblance to a Lam- 
bertian reflector. 


Experimental Procedure 

Prior to filling the tank with any water the bottom albedo, A^, was 
measured directly. The wavelength range considered for this measurement, as 
well as all others to be presented, was from 400nm to 700nm in increments of 
20nm. After measuring A^, several water types varying in absorptive capacity 
were added to the tank one at a time. The absorptance of each water type was 
controlled by adding known quantities of rhodamine dye. 

Table 1 is a summary of the physical characteristics associated with 
each water type considered. In each case, volume reflectance was calculated 
by normalizing the recorded upwelling radiance from the tank to that re- 
flected from the barium sulphate panel. Thus, 



0.79 L 

u 


ttL 


RP 


( 10 ) 
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where 0.79 is the reflectance of the barium sulphate panel. 


The only variable within equation (9) not directly measured was . 
Values were calculated by simultaneously solving equation (9) for two dif- 
ferent depths of each water type which yields 




( 11 ) 


Values chosen for the constant terms in equation (9) are as follows: 
= 0.06, = 0.02, Q = TT, and n^ = 1.33-1.39. The latter range of 

values reflects the dependence of upon the wavelength considered. 


Results and Conclusions 

Figure 3 is a plot of versus wavelength for each of the water types 

considered. With the addition of 2.7 ppm rhodamine dye to the tap water the 
value of increase sharply in the high absorption region between 480nm and 
580nm. A maximum value occurs at 560nm. From 580nm to 700nm, drops to 
near that of clear tap water. The same trend, only more exaggerated, occurs 
for tap water with the addition of 10.8ppm rhodamine dye; the width of the 
absorption band increases to include 460nm and 580nm, and the absorption peak 
occurs at 540nm rather than 560nm. 

In each water type with dye added, similar values of occur between 
600nm and 700nm. It is interesting to note that such values are lower than 
those representing clear tap water with no dye added. Clearly, the dye is 
transparent for longer wavelengths which accounts for its characteristic red 
tinge. Yet, intuitively, this window should not be any more transparent than 
the clear water to which the dye was added. If indeed equation (10) accur- 
ately describes K^, then the differing values in the longer wavelengths could 
be a result of the varying quality of the tap water used in the preparation 
of each water type. It was noted in the case of the clear tap water that the 
water used appeared to have a slight green tinge about it. As such, an 
associated absorption band in the red portion of the spectrum could account 
for the shape of the clear water curve. The tinge was not noticed in pre- 
paring either of the solutions containing rhodamine dye. 

According to equation (11), the total attenuation coefficient may be 
calculated in terms of the change in upwelling radiance with respect to the 
change in water depth. The calculated values of are thus independent of 
bottom albedo. Rewriting equation (9) in terms of the bottom albedo with 
Q = TT and ignoring the reflectance term yields the relationship 

ttL n^ d 

= E^(l-s")(l-p^y " (12) 
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Utilizing equation (12) , model predictions for were derived in terms of 
the calculated values of K^. As a comparison to the direct measurements of 
Aj, the calculated values of A^ as’ well as those directly measured were 
plotted with respect to wavelength and are presented in figure 4. 

Model predictions for A^ were found to be quite similar to the directly 
measured values. A statistical correlation coefficient, r, was calculated as 
a measure of the similarity between the two curves. The comparison of measured 
to predicted values of A^ resulted in r equivalent to 0.957, a significant 
correlation. 

The measured values of A, as well as the calculated values of were 
used in conjunction with equation (9) to calculate model predictions for the 
upwelling radiance from the tank. Several depths for each water type were 
considered. Figure 5 represents both predicted and measured values of 
radiance reflectance (above the water surface) from clear tap water for three 
different water depths. The predicted values were found to correlate well 
with the measured values. Correlation coefficients of 0.974, 0.964, and 
0.994 were calculated for the curves representing water depths of 3.81cm, 
22.86cm, and 76.52cm respectively. 

Comparable results were obtained for tap water with the additions of 
2.7ppm and 10.8ppm rhodamine dye and are shown in figures 6 and 7 respec- 
tively. Again, predicted values were quite close to the measured values. 

For the tap water with 2.7ppm rhodamine dye added, values for r ranged from 
0.963 to 0.983. In the case of tap water with the addition 10.8ppm rhodamine 
dye, values ranged from 0.959 to 0.979. 


QUALITATIVE PREDICTIONS 

We are now in a position to make some qualitative predictions concerning 
the way the reflectance of the water changes when the water is stratified. 

As an example, we will consider a situation which might well occur at a 
frontal boundary, such as at the outer edge of the Chesapeake Bay plume. 

The situation is illustrated in figure 8 which shows two adjacent water 
masses differing in color. A tongue of green water overlies a portion of the 
blue-green water. 

Each water mass was characterized optically using measurements published 
in reference 6. Station 4, in Apalachee Bay where the river effluent 
contains significant qualitites of dissolved organic material from the large 
inland swamp regions, was taken as characteristic of the green water. Station 
6b, in blue-green water near the mouth of the Mississippi, was chosen to 
represent the other water type. The irradiance reflectances and diffuse 
attenuation coefficients for these two stations are shown in figure 9. For 
this situation equation (7) becomes 

-K Az -K Az 

R = R- e ^ + R (1-e ^ ) (13) 

bg g ' 
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where R is the irradiance reflectance. The subscripts bg and g refer to 
blue-green (Sta. 6b) and green (Sta. 4) waters respectively. We have 
neglected the affects of bottom reflectance (d ->■ “) . Since, for the 
purpose of modeling, the layers were assumed to be of optically Infinite 
horizontal extent, the boundary between the two water masses is approximated 
in figure 8 by distinct levels. 

Results are shown in figure 10 where R(O-) is plotted vs. wavelength 
for several thicknesses of the upper layer. When Az=0 the reflectance 
is that of the blue-green water alone. As the thickness of the upper 
layer increases, the reflectance approaches that of optically deep green 
water. Note that at A=10m, the reflectance of the two-layer system is 
essentially indistinguishable from that of optically deep green water. 

At this depth in these waters, a highly reflective bottom would be 
easily visible; the blue-green water has little effect due to its low 
reflectivity. The effective penetration depth for these waters, for the 
purposes of remote sensing, is only ^lOm although a reflective bottom 
might be detectable at 2 or 3 times this depth. 

There are several points worth noting in comparing these curves. 


(1) The change in color is quite rapid; most occur for a layer, 
thickness of 2 attenuation lengths or less. This agrees with 
the conclusions of Gordon and McCluney (ref. 7). 

(2) The most obvious point about figure 10 is the existence of 
a nodal point at 480nm. At this particular wavelength the 
reflectance of both water types is the same 

even though the optical properties may differ significantly 
(see figure 9) . It is not possible to distinguish between 
the two water types by their reflectance at this nodal wave 
length. On the other hand, a nodal wave-length, when it 
exists, will provide an ideal point of reference when obser- 
ving a two-component, stratified system. 

(3) Less obvious, but at least as important, is the fact that the 
rate of change in reflectance is wavelength dependent. As 
can be seen from an examination of equation (13) the rate of 
change is entirely dependent on the spectral diffuse attenuation 
characteristics of the top layer which is the green water in 
this case. 


CONCLUDING REMARKS 

A relatively simple mathematical model of radiative transfer in a 
vertically inhomogeneous water mass has been presented. The model is 
quite simple in concept and is primarily designed to illustrate the ways 
in which water color may vary in situations which are difficult to model 
exactly. The preliminary experimental results presented above are 
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remarkably accurate, indicating that the model is faithful to reality in 
at least these simple situations. 

Also presented were some qualitative predictions relating to a 
situation which might well arise in highly dynamic estuarine regions 
such as the Chesapeake Bay Mouth. These predictions suggest that, when 
the water column is stratified, considerable variation in color might 
occur if the top layer is variable in thickness within two attenuation 
lengths of the surface. 
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Table 1 


Physical Parameters of Water Tank 


WATER TYPE 

TRIAL 

TIME 

WATER DEPTH 

BOTTOM TYPE 

False Bottom 

1 

1202 

0.00 cm 

Flat white 

Tap Water 

1 

1213 

3.81 

Flat white 


2 

1232 

22.86 



3 

1255 

49.85 



4 

1314 

76.53 


Tap Water 

1 

1322 

74.93 

Flat white 

With 2.7 ppm 

2 

1333 

50.83 


Rhodamine Dye 

3 

1343 

25.40 


Added 

4 

1350 

11.43 


Tap Water 

1 

1429 

76.20 

Flat white 

With 10.8 ppm 

2 

1442 

50.80 


Rhodamine Dye 

3 

1451 

26.04 


Added 

4 

1458 

12.70 
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Figure 3.- Plot of the two-way attenuation coeffi 
cient versus wavelength for each water type 
where is equivalent to the sum of 
equations (1) and (2) . 
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Figure 4.- Plot of direct measurements and model 
predictions of the reflectance of the 
tank bottom versus wavelength. 



WAVELENGTH (nm) 

Figure 5.- Plot of measured and predicted values 
of voluTiie reflectance versus wavelength for the 
case of clear tap water - 



WAVELENGTH (nm) 

Figure 6.- Plot of measured and predicted values of 
volume reflectance versus wavelength for the 
case of clear water +2.7 ppm Rhodamine dye. 
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WAVELENGTH (nm) 

Figure 7.- Plot of measured and predicted values of 
volume reflectance versus wavelength for the 
case of clear water + 10. .8 ppm Rhodamine' dye. 



Figure 8.- Boundary between two water types. Green 
water overlies blue-green water. The depth is 
scaled to the diffuse attenuation length for 
the green water at 550 nm. 






